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Abstract 

We report on the progress in work on improving precision of the standard model theoretical 
predictions for the top quark pair production and decay into six fermions at a linear collider. Two 
programs have been combined into a single Monte Carlo program: eett6f , a MC program for 
e + e~ — > 6f, and topf it, a program for electroweak radiative corrections to e + e~ — * ti. The MC 
program is described and preliminary numerical results are shown. 

PACS numbers: 12.15.-y, 13.40.Ks, 14.65.Ha 
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I. INTRODUCTION 



In order to disentangle possible new physics effects that may be revealed in the process 
of top quark pair production 

e + e" -> tt (1) 

at a linear collider from predictions of the standard model (SM), one needs to know the 
latter with high precision, possibly at the level of a few per mill [lj]. That high precision 
level requires taking into account radiative corrections. Because the top (and antitop) quark 
of reaction (0) decays into a b (and b) and fermion-antifermion pairs, we should study the 6 
fermion reactions of the form 

e + e~ - bhKbhft, (2) 

where fi, f 2 = u e , v r , u, c and f[,fz = e~, t~ , d, s. 

Reactions (J2J) receive contributions from typically several hundred Feynman diagrams 
already at the tree level, e.g.the hadronic reaction e + e~ — ► budbdu in the unitary gauge, 
neglecting the Higgs boson couplings to fermions lighter than a 6-quark, gets contributions 
from 1484 Feynman diagrams. This is the main reason why the calculation of the full 0(a) 
radiative corrections to any of reactions (J2J) seems not to be feasible at the moment. 

Therefore, in order to improve the precision of the lowest order predictions one should 
try to include at least the leading radiative effects, such as initial state radiation (ISR) and 
factorizable radiative corrections to the process of the on-shell top quark pair production 
(JTJ) and to the decay of the top and antitop quarks. 

In the present lecture, we report on the progress of the SM theoretical predictions for the 
top quark pair production and decay into six fermions at a linear collider, paying special 
attention to the precision improvement. 

We will show a sample of results on reactions (121) including some higher order effects 
including the ISR correction obtained with eett6f |2J, a Monte Carlo program for reactions 
e + e~ — > 6 fermions relevant for the top production. Further, we introduce the one-loop 
electroweak corrections to reaction (JTJ obtained with topf it a numerical package for 
calculating electroweak radiative corrections to the on-shell pair production of a massive 
fermion-antifermion pair. The two programs have been combined into a single Monte Carlo 
program. 
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II. CALCULATING ISR WITH EETT6F 

The ISR correction has been implemented recently into eett6f 4] in the leading loga- 
rithmic (LL) approximation utilizing the structure function approach. In this approach, the 
corrected differential cross section da LL (pi,p 2 ) of any reaction (J2J) reads 

do-Bom+isR {pi,P2) = J dxi d x 2 (x x , Q 2 ) (x 2 , Q 2 ) 

> : -da Born (x 1 p 1 ,x 2 p 2 ) , (3) 

where pi (p 2 ) is the four momentum of a positron (electron), x\ (x 2 ) is the fraction of 
the initial momentum of the positron (electron) that remains after emission of a collinear 
photon and dt7Bom(xipi,x 2 p 2 ) is the lowest order cross section calculated at the reduced 
four momenta of the positron and electron. The structure function (x, Q 2 ) is given by 
Eq. (67) of [5], with 'BETA' choice for non-leading terms. The splitting scale Q 2 , which is 
not fixed in the LL approximation is chosen to be equal s = (pi + p 2 ) 2 . 

with the 
where 

also the corresponding lowest order total cross sections 0"B O m are compared. The reader 
is referred to |2j for a comparison of the lowest order predictions with other existing MC 
programs. The input parameters are the same as in [f| and cuts are given by 

9(1, e ± ) > 5°, 9(q,e ± ) > 5°, 9(1,1') > 5°, 9(1, q) > 5°, 

Ei > 10 GeV, E q > 10 GeV, m(q, q') > 10 GeV, (4) 

where 9(i,j) is the angle between particles i and j in the laboratory system, q and / denote 
a quark and a final state charged lepton, respectively, and m(q, q') is the invariant mass of 
a qq' quark pair. 

The results of both programs agree nicely, basically within one standard deviation, except 
for e + e" — > bv e e + duh, where the deviation is bigger, amounting to about 3 standard 
deviations. The discrepancy should be clarified by fine tuned comparisons. Differential 
cross sections obtained with eett6f that are plotted in FigG] agree with those presented in 
within accuracy of the plots. 



The ISR corrected cross sections, cxBorn+isR , of different channels of (J2J obtainec 
current version of eett6f are compared with the results of Lusiper Q in Table 1 ' 
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TABLE I: Comparison of the cross sections of different channels of (J2J) at yfs = 500 GeV of Lusifer 
[(| and eett6f 3, J]. The input parameters and cuts, see Eqs. @ are the same as in 6]. All cross 
sections are in fb. The number in parenthesis show the uncertainty of the last decimals. 



_l_ 

e^e — > 


LUSIFER 


eett6f 


°"Born 


0"Born+ISR 


CBorn 


CBorn+ISR 


bv e e + e~u e b 


5.8530(68) 


5.6465(70) 


5.8622(63) 


5.6441(67) 


bh' e e + fi^u fl b 


5.8188(45) 


5.6042(38) 


5.8189(37) 


5.6075(59) 


bv^^ iTv^b 


5.8091(49) 


5.5887(36) 


5.8065(33) 


5.5929(54) 


b [i + t~9 t b 


5.7998(36) 


5.5840(40) 


5.7992(32) 


5.5844(33) 


bv^ fi + dub 


17.171(24) 


16.561(24) 


17.213(23) 


16.569(17) 


without QCD: 


17.095(11) 


16.454(10) 


17.106(15) 


16.457(16) 


bu e e + dub 


17.276(45) 


16.577(21) 


17.301(26) 


16.741(43) 


without QCD: 


17.187(21) 


16.511(12) 


17.149(16) 


16.522(17) 



III. CALCULATING ONE-LOOP ELECTROWEAK CORRECTIONS WITH 

TOPFIT 

topf it 0] is a numerical package for calculating the complete one-loop corrections to 
top-pair production JI| at a linear e + e~ collider in the continuum energy region. The 
corrections are represented in terms of six independent form factors, which is suitable for 
implementation in the Monte Carlo generator. We recall the normalization of the differential 
cross section 



da(e + e-^ti) f3 t 
d cosd 327r s 



[\M B \ 2 + 2$te{M* B 5M) 



(5) 



conf 



where $ is the antitop production angle with respect to the initial e + beam, (3 t = 

1 /2 

(1 — 4 mf I s) and the sum stands for the spin and color configuration of external fermions. 
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FIG. 1: Invariant mass (left) and angular (right) distributions of a Mu-quark triple at y/s = 500 
GeV with cuts. 



There are four nonvanishing form factors Ffg in Born approximation 

M B = E Ff B M hab , (6) 

o,6={l,5} 

while the one-loop correction 5M. to the amplitude can be expressed in terms of six inde- 
pendent form factors: F± b , a,b = {1,5}, plus additionally F 3 X1 and F 3 51 : 

SM = £ Ff M ltab + F^ M^n + Fl 1 M 3M . (7) 

a,6={l,5} 

The invariant amplitudes in Eqs. (jOJ) and (J2J) are given by 

iMi,ab = [v(pi) Y9au(p 2 ) ] [u(p 3 ) l^gbV^Pi) ], a, 6 = {1,5}, 

^3,11 = ~[v(Pi) ^s1«(P2) ] [«(pa) lv(P4) ], (8) 

^3,51 = -[«(Pl) ^3 75^(^2) ] [«(P3) 1«(P4) ]■ 

The normalization of the form factors is chosen such that the pure photonic Born ampli- 
tude becomes F^ 1 = e 2 Q e Q t /s; (see j^] for the details on the analytic expressions). 

The resulting corrected differential cross section formula can then be written in terms of 
the computed form factors, the top quark color factor (c t ) and the kinematical Mandelstam 
variables s,t and u (s + t + u = 2m 2 ): 
da 2 



d COS 7? s s 



f^^eUsm 2 (F?F$ - F?F% + F?F% - F?F% 



TABLE II: Born and weak one-loop form factors for e + e — > tt at ^/i = 500 GeV and cost? = 0.7. 





Born F.F. (GeV~ 2 ) 
3?e 9m 


Weak one-loop F.F. (GeV~ 2 ) 
3?e 9m 


psi 

Ff 
m t Fl l 
mtF^ 1 


-2.50926472- 10~ 7 
1.56200966- 10" 8 
5.62400131- 10~ 8 
-1.37479846- 10~ 7 




1.08874859- 10~ 8 -2.55015622- 10~ 9 
-9.75979648- 10~ 9 -9.02716758- 10~ 9 
-7.00614109- 10~ 9 -6.39336661- 10~ 9 
-1.11566191- 10 -9 7.86314816- 10" 9 

8.14431322- 10~ 10 -8.44590997- 10~ 10 
-9.09997715- 10~ 10 4.88604910- 10~ 10 



_ 2 ( M 2 - t 2 - 2m 2 (u - t)) (F^F^ + FfFg + F^F^ + F^F, 



51» 
IB 



+2{u 2 + t 2 + 2m 2 (s - m 2 )) ]T F^Ffg 

a,6={l,5} 



+A(ut - rat) {Fs'Fib + ^ib ) }, (9) 
where the dimensionless barred form factors are defined for convenience as 



— s 
*1B = ~2*1B > 



Ff ee ±(i F ^ + Ff), (10) 
s 

7 2 



Ff = -m t Ff. 



topf it is rather flexible and includes also real photonic bremsstrahlung. In order to fit 
the needs of the MC program described here, some specific values of flags have been chosen: 
I WEAK = 1 (inclusion of pure weak corrections), IQED = IQEDAA = 0. The latter two settings 
switch off the calculation of photonic corrections due to bremsstrahlung (IQED, with the 
related virtual corrections) and of the fermionic vacuum polarization effects (IQEDAA). The 
sample corrections were computed using the same input parameters as in ^]. In Table ILT1 
form factors are given for y/s = 500 GeV and cos$ = 0.7, and in Table ITTT1 some differential 
cross section values for ^fs = 500 GeV. 
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TABLE III: Differential cross sections for e + e — ► tt at = 500GeV. 



COS 17 


do- 


Born (pb) 


do- 


Born+Weak lP ' 


d cos 1? 


dcos 1? 


-0.9 


0.1088391941 




0.1011751626 


-0.7 


0.1218770828 




0.1136440560 


-0.5 


0.1422750694 




0.1332668988 


0.0 


0.2254704640 




0.2115640646 


0.5 


0.3546664703 




0.3300633143 


0.7 


0.4192250441 




0.3883965473 


0.9 


0.4911437158 




0.4528185418 



IV. NUMERICAL RESULTS 

In the present section, we will show the results of implementation of the one-loop 
electroweak form factors for on-shell top quark pair production obtained with topf it in 
the MC program eett6f. Amplitudes (JHJ) have been implemented in eett6f with the 
helicity amplitude method of 0, EiJ- The actual numerical values of the form factors 
Fib,Fi , a,b = {1,5}, -F 3 n and _F 3 51 are computed by topf it and then transferred to 
eett6f . Unfortunately, the direct calculation of the one-loop corrections slows down the 
MC computation by more than a factor 1000. In practice, this problem has been solved 
applying the following strategy: the values of the one-loop form factors are computed for a 
given c.m.s. energy for a grid of several hundred, or a few thousand different values of cosine 
of the antitop scattering angle, cosi9. These fixed values of the form factors are stored in the 
computer memory and then used for computation of the form factors at any intermediate 
value of cos d by means of a linear interpolation between the neighbouring fixed values. The 
interpolation routine has been tested and shown to be very precise and efficient. 

As a first step we calculate the differential cross section of reaction (J2J) in the narrow 
width approximation for the top and antitop, which is obtained by multiplying the one-loop 
corrected cross section (JHJ) with the corresponding branching fractions leading to a specific 
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final state of (J2J) 

da (e+e- - £t -> bfJlbfJl) =da(e + e - -> tt) 

x dr (t - 6/x/jf) dr (t^ 6/2/2) A?- (11) 

This approach is gauge invariant and can be used as the reference for any other way of 
implementation of the one-loop electroweak corrections in the MC generator for a simulation 
of reactions (J2J) to (9(a). 

Typical results for one specific channel of reaction (J2J in the narrow top width approx- 
imation of Eq. are shown in Figs. El and H3 The differential cross sections of reaction 
e + e _ — > bv^n+bdu at a/s = 500 GeV without cuts are plotted on the left hand side of Fig. El 
Both the lowest order (dotted line) and one-loop corrected (solid line) differential cross 
sections in the 6-quark energy are shown. The same differential cross sections with cuts of 
conditions (jlj) are plotted on the right hand side of Fig. El The corresponding differential 
cross sections in the energy of /x + are plotted in Fig. El 



e + e — > tt — > bu^fi^bdu 
Born 

Born + weak 



y/a = 500 GeV 



An 
AE h 



100 150 
E b (GeV) 



200 250 



£ 


1 1 1 1 

+e~ -J- ti bv^fi+bdu ys = 500 GeV 




Born 








(with cuts) 















50 100 150 

E b (GeV) 



200 250 



FIG. 2: Energy distributions of a 6-quark at y/s = 500 GeV without cuts and with cuts. 

We see that the quite substantial effect of the one-loop electroweak correction on differ- 
ential cross sections of e + e~ — > bu^fi + bdu at y/s = 500 GeV is reduced by the cuts (JH). This 
kind of reduction is of course not desirable if one wants to test the consistency of the SM 
at the quantum level. However, it might be helpful in performing relatively fast numerical 
simulations of different channels of reaction (J2J with the use of the MC generators working 
only with the lowest order amplitudes. 
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FIG. 3: Energy distributions of a ji + at y/s = 500 GeV without cuts and with cuts. 



V. SUMMARY AND OUTLOOK 



We have reported on the progress in our work on improving precision of SM theoretical 
predictions for the top quark pair production and decay into six fermions at a linear collider. 
Apart from including the ISR effects in eett6f |4J, we have created a MC program that 
combines two programs eett6f and topf it. The program allows for including one-loop 
electroweak corrections to the top quark pair production in reactions (J2J in the narrow 
width approximation for the top and antitop. The description of these corrections to the 
on-shell top pair production has to be complemented by the corresponding ones to the 
decays of the too q uarks Q, in otdet to complete the doub,e pole approbation. This is 
in preparation [12J. 
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